;Geometric double-quantum singlet excitation sequence.
;
;$CLASS=HighRes
;$DIM=1D
;$TYPE=
;$SUBTYPE=
;$COMMENT=
;$RECOMMEND=y

#include <Avance.incl>
#include <Delay.incl>
#include <Grad.incl>

"acqt0=-p1*2/PI"
"p2=2*p1"
"d15=2*d7"
"p4=0.5*p1"
"d4=d7"

1 ze 
  

2 30m


;SOD filter

;T00 filter 
 10u LOCKH_ON
3 5u pl1:f1 
  40u UNBLKGRAD 
  p16:gp1 
  2m  
  (p1 ph4):f1 
  p19:gp2 
  2m 
  (p1 ph4):f1 
  (p1 ph5):f1 
  p29:gp3
40u BLKGRAD
  2m 
;JSE train of echoes
4	d4
	p1 ph1
	p2 ph1+90
	p1 ph1
	d4
lo to 4 times l20
lo to 3 times l2
  10u LOCKH_OFF

5 d1                       


;First GEODQ block

6	p1 ph2

7	d15
	p1 ph1
	p2 ph2
	p1 ph1
lo to 7 times l1

8	d7
	p1 ph1
	p2 ph2
	p1 ph1
	d7
lo to 8 times l1

9	p1 ph2		
					
 
;S2DQ

10	p1 ph7                

11	p4 ph6                 

12	d4                     
	p1 ph6
	p2 ph2
	p1 ph6
	d4
lo to 12 times l1

 
 ;T00 filter
  5u pl1:f1 
  40u UNBLKGRAD 
  p16:gp1 
  2m  
  (p1 ph4):f1 
  p19:gp2 
  2m 
  (p1 ph4):f1 
  (p1 ph5):f1 
  p29:gp3
  40u BLKGRAD
  2m 

    
     2m						            	 ;Evolution delay
 

;T00 filter
  5u pl1:f1 
  40u UNBLKGRAD 
  p16:gp1 
  2m  
  (p1 ph4):f1 
  p19:gp2 
  2m 
  (p1 ph4):f1 
  (p1 ph5):f1 
  p29:gp3
40u BLKGRAD
  2m 

;S2DQ

14 	d4                     
	p1 ph6
	p2 ph2
 	p1 ph6
	d4
lo to 14 times l1

15	p4 ph6            

16	p1 ph7            

;Second GEODQ block

17	p1 ph20

18  	d15
    	p1 ph3
	p2 ph20
	p1 ph3
lo to 18 times l1

19	d7
	p1 ph3
	p2 ph20
	p1 ph3
	d7
lo to 19 times l1


  go=2 ph31 

  30m mc #0 to 2 F0(zd)
  lo to 1 times td1

exit


;Phases
;------T00 filter phases--------
ph4 = (65536) 9965 9965  
ph5 = 2 2 
;------Pulse phases--------
ph1 = 0 0 0 0 		
ph2 = 1 														
ph6 = 0 														
ph7 = 3 														
ph3 = 0 1 2 3 		
ph20 = 1 2 3 0 											
;------Receiver phase--------
ph31 = 0 3 2 1 


;pl1 : f1 channel - power level for pulse (default)
;p1 : f1 channel -  90 degree high power pulse
;p2 : 2 * p1; f1 channel - 180 degree high power pulse
;d1 : relaxation delay; 1-5 * T1
;d11: delay for disk I/O                             [30 msec]
;ns: 1 * n, total number of scans: NS * TD0
;l1: Echo number: PI/(2*arctan(D/J))
;l20: Singlet Order Destruction Echo Number
;l2: SOD cycles
;d15: Delay: PI/omega
;d7: Delay: PI/(2*omega)

;omega: 2*PI*SQRT(J^2+D^2)
